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Abstract 



A system of periodic poly-disperse coated nano-rods is considered. Both the coated nano- 
rods and host material are non-magnetic. The exterior nano-coating has a frequency dependent 
dielectric constant and the rod has a high dielectric constant. A negative effective magnetic 
permeability is generated near the Mie resonances of the rods while the coating generates a neg- 
ative permittivity through a field resonance controlled by the plasma frequency of the coating 
and the geometry of the crystal. The explicit band structure for the system is calculated in 
the sub-wavelength limit. Tunable pass bands exhibiting negative group velocity are generated 
and correspond to simultaneously negative effective dielectric permittivity and magnetic perme- 
ability. These can be explicitly controlled by adjusting the distance between rods, the coating 
£h ■ thickness, and rod diameters. 

Short title: Tunable Double Negative Band Structure 
PACS 42,78 

> 
00 ; 1 Introduction 

Metamaterials are new class of man made materials that impart unconventional electromagnetic 
properties derived from sub- wavelength configurations of different conventional materials pQ. The 
l/"~j ■ first such materials were seen to exhibit behavior associated with negative bulk dielectric constant [5] 

and were constructed from a cubic lattice of metal wires. Subsequently negative effective magnetic 
permeability at microwave frequencies [T] were derived from periodic arrays of non-magnetic metallic 
split ring resonators [1 . Double negative or left handed metamaterials with simultaneous negative 
bulk permeability and permittivity at microwave frequencies have been developed using periodic 
arrays of metallic posts and split ring resonators [3]. Subsequent work has delivered several new 
designs using different configurations of metallic resonators for double negative behavior 4, 5, 6, 7 '. 

ei ia ma . 

For higher frequencies in the infrared and optical range an alternate strategy for constructing 
negative effective permeability from non-magnetic components relies on Mie resonances associated 
with small rods or particles made from high permittivity materials |11) . A double negative meta- 
material may be achieved by coating the high permittivity material with a frequency dependent 
plasmonic or Drude type coating |12j as well as coatings with more general frequency dependence 

In this article we focus on the second approach and propose periodic assemblages of aligned 
non-magnetic dielectric rods each clad with a non-magnetic plasmonic coating. For this case we are 
able to explicitly calculate the propagation band structure in the sub-wavelength limit. The pass 
bands and stop bands are given by formulas that depend explicitly on the rod diameters and coating 
thickness. We show how to construct tunable pass bands associated with double negative effective 
properties. The physical origin of the negative effective permeability is due to the excitations of 
Mie resonances inside the rods. On the other hand the negative effective permittivity is caused by 



the extreme dielectric properties near the plasma resonance of the coating. We provide the explicit 
relationship linking Mie resonances and the frequency dependent effective permittivity to the spacing 
of the rods, the rod radii, and the coating thickness for the coated rod assemblage proposed here. 

These relationships are found not by appealing to effective medium theory based on the Clausius- 
Mossotti formula, but instead we show that it is profitable to take a different approach and charac- 
terize wave propagation using explicit multiscale expansions for waves inside the photonic crystal. 
The wave number associated with a Bloch wave inside the d-periodic crystal is denoted by k — 2ir/\ 
where A is the wavelength. The approach taken here provides an explicit power series expansion for 
the fields in the parameter r\ = dk. We outline a systematic framework in which the homogenized 
dispersion relation is recovered directly from the expansion in the subwavelength d « A limit. 
We introduce photonic crystals made from coated rod assemblages and obtain explicit homogenized 
dispersion relations for these geometries. The dispersion relations provide explicit conditions on the 
the distance between neighboring rods, rod radii, and coating thickness necessary for generating pass 
bands associated with double negative behavior. Two examples are provided showing how the band 
structure can be manipulated to generate double negative pass bands with negative group velocity. 

We conclude referring to related work addressing wave propagation inside high contrast media. 
When the crystal period and wavelength have the same length scales, the ability to open band gaps 
for photonic crystals is established and developed in the set of papers [14] , [15] , and [16] . The power 
series strategy presented here for subwavelength analysis has been utilized earlier and developed in 
[17j for characterizing the dynamic dispersion relations for Bloch waves inside plasmonic crystals and 
for understanding the influence of effective negative permeability on the propagation of Bloch waves 
inside high contrast dielectrics |18j . Earlier work in two scale homogenization theory for high contrast 
dielectrics delivers a frequency dependent effective magnetic permiability [19], [20] identical to the 
generic effective permeability tensor recovered here, see also the work of [21 for homogenization 
in high contrast media. The connection between high contrast interfaces, homogenization and the 
associated generation of negative effective magnetic permeability is made in [22j . More recently 
two-scale homoginization theory has been developed for three dimensional split ring structures that 
deliver negative effective magnetic permeability |23] and for metal fibers delivering negative effective 
dielectric constant [24] . A method for creating metamaterials with prescribed effective dielectric 
permittivity and effective magnetic permeability at a fixed frequency is developed in 25 . 

2 Electromagnetic fields inside photonic crystals made from 
coated rods and the subwavelength limit 

The meta material is a two-dimensional photonic crystal made of parallel coated rods, see Figure [T] 
There can be one or more coated rods inside the crystal period. The time harmonic field is p- 
polarized and the magnetic field inside the crystal is H = H(x) exp (— iwt)e3 where x = (£1,0:2) in 
the a;ia;2-plane. 




Figure 1: Photonic crystal made from coated rods 



The period for the crystal is d and the dielectric coefficient e<z(x) takes the values 



£d( x ) 



(R 


in the rod, 


e p 


in the coating, 


th 


in the host material. 



(2.1) 
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The coating is a cylindrical shell of plasmonic material with dielectric constant e p (w 2 /c 2 ) = 1 — Jjw g ■ 
Here uj p is the plasma frequency associated with the coating material and c is the speed of light in 
vacuum. The dielectric constant of the rod is chosen according to 6r — ^h where e r has units of area. 
The idea is to choose the dielectric permittivity of the rod to be large so that the corresponding Mie 
resonances are excited in the sub- wavelength limit. The dielectric constant of the host material is 
given by e h = l. 

The time harmonic magnetic field for the d-periodic crystal is a Bloch wave 



H = /i(x) exp (ikk ■ x), 



(2.2) 



where h is d-periodic in x, the wave number is given by k, and the direction of propagation in the 
xiX2-plane is described by the unit vector k = («i, k 2 ). H satisfies the Helmholtz equation 



vVv# = - -v#, 



(2.3) 



with propagation bands described by dispersion relations involving ui and kk. We examine the 
situation in the sub- wave length limit when d tends to 0. Instead of developing an effective medium 
theory based on the Clausius-Mossotti formula we show that it is profitable to take a different 
approach and expand the Bloch wave in a multiscale power series in r\ = kd. Here we develop a 
multiscale power series expansion for h to recover "homogenized sub- wavelength" dispersion relations 
for plane waves inside a magnetically active double negative effective medium. 

We introduce the variable y = x/d that takes the period cell of size d into the unit period cell 
Y. Here Y is split into the subdomains given by the union of rod cross sections R, the rod coatings 
P and the host H. Figure illustrates a unit cell Y containing one coated rod cross section. The 
expansion parameter is given by rj — kd = 2 ^r L - The sub- wavelength limit is described by r\ — > for 



A fixed. In what follows it is useful to introduce £ defined by £ 2 
satisfied by h for y in the unit period Y are 



c 2 k 2 



Lh 
-Lh 



rj 1 ^ 2 h 1 in the host, 
r] 2 ^ 2 ep(^ 2 k 2 )h, in the coating, 
e r £, 2 k 2 h, in the rod, 



The propagation equations 



(2.4) 



where 



Lh 


= Ah + 2i-qk ■ Vh - ifh, 


e P (ek 2 ) 


u 2 P /c 2 

ek 2 ' 


the transmission conditions 




n ■ (V/i + irjkh) 


= n-e p 1 (£ 2 fc 2 )(V/i + i?7K/ 



l {( 1 k 2 ){Vh + ir 1 kh) 



k 2 ( 



a 



(V/i + irjkh) 



, H-C interface, 
R-C interface. 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



Here the interface separating the host phase and coating phase is denoted by H-C and the interface 
between the coating phase and the rod is denoted by R-C. Evaluation of quantities on the host side 
of the interface is denoted by the subscript H , on the coating side by C, and by R on the rod side. 
The unit normal vector n on the host - coating interface points from the coating into the host and 



the unit normal vector on the rod - coating interface points from the rod into the coating. We 
expand h and £ 2 as 



h(dy) = ho(y) + r)hi(y) + v 2 h 2 (y) 

e = eo+vei+v 2 e 2 + --- 



(2.9) 
(2.10) 



where each term in the expansion is continuous in y and periodic on Y . Substitution of the expansion 
into (|2.4[) , (J2.7P , and (|2.8[) equating like powers of r] gives the following equations used to determine 
ho- 



— Aha = 0, in the host, 
—Aha — 0, in the coating, 
— Aha = e r£ok h, in the rod, 



and transmission conditions, 

n ■ (VV 



H 



ep 1 ^ 2 ) (V/io) , H-C interface, 



n-e p 1 (^k 2 )(Vho) = 0, R-C interface. 



(2.11) 
(2.12) 
(2.13) 



(2.14) 
(2.15) 



Noting further that h is continuous across these interfaces we apply (|2.11|) , (|2.12|l , (|2.14j) , and (|2.15|) 
to discover that ha is a constant function ho = h outside the rods. From linearity it follows that the 
ho field inside the rod can be written as ha — hm(y) where 



Am = er£,a^ m -> m the rod, 



(2.16) 



and m = 1 on the boundary of the rod. Define the function g(y) to be 1 for y outside the rod and 
g(y) to be m(y) for y inside the rod then ho(y) is defined up to a multiplicative constant h and is 
given by 



ho(y) = My)- 



(2.17) 



.'72 




2/1. 

Figure 2: Unit cell containing a single rod cross section. R represents the rod cross section , P the 
plasmonic coating and H the host. 



Next we determine h\ exterior to the rods. On equating like powers of r\ we find that 

-Ah x = 0, (2.18) 



outside the rods, and the corresponding transmission conditions for hi are given by 

n-(\/hi+ikh) = n ■ ep 1 (^k 2 ) (V 'hi + ikh) , H-C interface, (2.19) 

Ih c 

n ■ ep 1 ^ 2 ) (V/ii + ikh) = 0, R-C interface. (2.20) 

\c 

The equation (J2.18I) . transmission condition Q2.19P and boundary condition (|2.20[) determine h\ up 
to an additive constant C\ and from the linearity of the problem we write h\ = ihW ■ k + C\ , Here 
W(y) = (wi(y), W2(y)) where Wi, i = 1,2 are the periodic solutions of the cell problem 

— Awi = 0, for y outside the rod. (2-21) 

n-iywi + e) = n-e p 1 (Cofc 2 )(Vu; 4 + e l ) , H-C interface, (2.22) 

Iff " \c 

n-e p 1 (^ 2 k 2 )(\7w l + e l ) = 0, R-C interface, (2.23) 

c 

where e 1 = (1, 0) and e 2 = (0, 1). We define the region inside the period cell exterior to the rods by 
D and introduce the position dependent dielectric constant e(y) taking the value 1 inside the host 
and ep(k 2 t; 2 ) inside the coating. For future reference we introduce the tensor e~} f (£o^ 2 ) defined by 



, |/7 , M1 ,- J '- .-„wv-TT- r--/„ 



rtAfa 2 ) = / £ (y) ^W + I) dy, (2.24) 



D 



where / is the 2x2 identity tensor. 

Now we identify the effective magnetic permeability. Following pQ, |llj the average component 
of the B field parallel to the rods is given by 

B ave = -]L f h(x.)e l ^ x Ux 1 dx 2 . (2.25) 

" JY d 

Here Yd is a period cell for the d-periodic crystal. Substituting the expansion (J2.10J) for h into (|2.25[) , 
retaining the lowest order terms and taking d — > gives 

lim Bave = h g(y)dyidy 2 - (2.26) 

d^o J Y 

The average component of the H field parallel to the rods is given by [1] , [11] 

i f (0,0,d) 

H ave = - / h(*)e< kk -^d Xl dx 2 , (2.27) 



d ^(o,o,o) 
here the average is taken over any line parallel and exterior to the rods. Sending d — > gives 

lim H am = h. (2.28) 

d->0 

Using (J2.26P and (|2.28|) the effective permittivity is given by 

l im |^i =Me//(? 2^). (2.29) 

We recover the explicit form of the effective permeability A 1 e//(Co^ 2 ) f° r a distribution of circular rods 
inside the unit period cell. To do so recall the following Dirichlet eigenvalues E and eigenfunctions 
ip for each rod cross section given by 

— Atp = Eij) for y inside rod and ?/> = for y on rod boundary . (2.30) 

For rods with circular cross section of radius < R r < 1 and polar coordinates < p < R r , 
< 9 < 2-7T, let Xnm be the mth zero of the nth Bessel function, then the Dirichlet eigenvalues are 



Enm = {Xnm/Rr) 2 , and the corresponding eigenfunctions are { J n ( 2 ^ lrL p)e me , J„(^ L p)e -4 ™ 9 }. The 
associated normalized eigenfunctions are denoted by ip^; m . Notice that the eigenvalue Eo m has only 
one associated normalized eigenfunction i/; 0m = pl° l "!,f / "',, and (V'om} 2 = ^p- when n = 0. Here 

(•} denotes the area integral over a rod cross section. It is obvious that {4>nm) 2 = f° r n 7^ 0. 
Expanding m(y) in the basis VWi noting that to — 1 = outside the rods gives for each rod 



c 2 — 
For a distribution of rods with radii R r j, j = !,...,£, we have 



1 - E g2 i _ j^ onOVw (2.31) 



/• °° 47ri? 2 £ 2 P 

/M/C^ 2 ) = / 9(y)dy = 1 - £ JV(^;) £ V 2 r %2,i°„p , (2-32) 

Jy i=0 m=0 Xom ?0^ K jm 

2 

where N(R r j) is the number of rods with radius R r j and A;^ m = x "^ . 

Last we find the homogenized sub- wavelength dispersion relation. Recall that £ 2 fc 2 = ^ and 
from (l2~T0l) 



— = -J+V~j+V "f + -". 2 - 33 



where 



w 



2 



c- 



2 



= Zfkf, fori = 0,1,.... (2.34) 



The sub- wavelength dispersion relation identified as the dispersion relation between ujq and kk given 

2 
by ^§- — ^ 2 k 2 . We now recover this dispersion relation by deriving the formula for £q. Equating like 

powers of r\ gives the following equations for hi exterior to the rods. 

- (Ah 2 + 2ir]k ■ Whi - ifh) = ffi, in the host, (2.35) 

-(Ah 2 +2iT]k-Vhi-rj 2 h) = $ e P (£(, k 2 )h, in the coating, (2.36) 

and the transmission conditions 

( ilk 2 - *% ] n . (y h2 + ikhl ) + (^ 2 fc 2 ) n . (y hl + j£/Tj 

= (£o fc2 ) " ' (V/i2 + i«/n) + (C?^ 2 ) n ■ (V/ii + ikh) , H-C interface, (2.37) 

O lc 



(£ofc 2 )n- (V/12 +ikhx) + (£, 2 k 2 )n- (Vhx+ikh) + (gk 2 ) n -Vh 



lc 



= 7^7— I Cn fc2 ~ ^% I n • Vfto , R-C interface. (2.38) 

k z e r \ c l j <b 

Integrating and adding (12.35[) and f|2.36[) gives the solvability condition 

£0 / Tidy = - / (A/ia + 2»k • V/n - h) rfy 

- / e_1 (Cofc 2 ) (A/i 2 + 2i« • V/11 - 7i) dy. (2.39) 

On integrating by parts and applying (|2.13[) together with the transmission conditions (|2.15[) . (|2.19p . 
(|2~2"0"j) . (|2~37|l . and (j2~35f shows that (|2~39"|l is equivalent to 

Co = l*7f,($k 2 )e$f($k*)K • «. (2.40) 



Writing w 2 /c 2 = £ 2 fc 2 and substitution into (I2.40J) delivers the homogenized sub-wavelength disper- 
sion relation 



Wn 



J=n^, (2.41) 

where the effective index of diffraction n 2 , f is given by 

«e//=Me//(^)(e e 7 / (^)«-«) \ (2.42) 

Applying (12.91) and collecting results shows that to lowest order the Bloch waves H = /i(x)e^ feKX ~ tw ' 
inside the photonic crystal are given by the sub-wavelength expansion 

H(x,t) = (%(^) + VM^) + 0{v 2 )) expi(kk ■ x - t(u + rjUx + 0(r? 2 )). (2.43) 

In the subwavelength limit the spatial averages of H correspond to spatial averages of plane 
waves associated with a magnetically active effective medium. To see this take any planar averaging 
domain S and pass to the d — > limit in the average to get 

lim — --— / H(x,t)dxidx2 = lim- —-- / hg(—)exr>i(kk-x — tuJo)dxidx2 

d-*o Area(b) J s d^o Area(b) J s d 

1 f 

= — T-77- / /i e ffhexpi(kk ■ x — tuJo)dxidx2- (2.44) 

Areayo) Jg 

In a similar way, taking averages over any interval a < x% < b on any line parallel to the X3 axis not 
intersecting the coated rods gives 

"(0,0,6) 1 /-(0,0,&) 



lim- / H(x,t)dx3 = lim- / hg(— )expi(kk • x — tu>o)dx3 

d^Q b-a 7(o.o, a) d ^° o-a J (0,0,0) « 

/■(0,0,6)_ 

■ / hexpi(kk ■ x — tuio)dx3. (2.45) 

J(0,0,a) 



b — a 



Thus the appropriate averages of the plane waves Hh om — hexpi(kk ■ x — tuj$) and 

Bhom = Heffhexpi(kk ■ x — two) provide approximations to the average field seen in the d periodic 

photonic crystal for < d << k. 

3 Electromagnetic fields inside coated rod assemblages 

Here we introduce a special class of photonic crystals made from coated rod assemblages and derive 
an explicit formula for the subwavelength dispersion relation (|2.4ip . The formula shows how the 
band structure depends explicitly on the distribution of rod radii, coating thickness and distance 
between neighboring rods. The assemblage is described as follows. Consider first the unit period cell 
Y filled with disks with radii ranging down to the infinitesimal. Inside each disk we place a centered 
rod cross section then a concentric coating of plasmonic material and finally a concentric coating of 
host material, see Figure El 

The cross sections of the coated rods together with the connected host material is depicted in 
Figure |U The ratio of the radii of the rod cross section, plasmonic coating, and host coating are the 
same for all the disks and we denote the area fractions of the host, coating, and core phases by 9h, 
dp , and Or respectively, see Figure [3] The radii of the rods used in the assemblage is denoted by 
R r \ > R r 2 > R r 3 • ■ ■ • The number of rods having radius R r j are denoted by N(R r j). For rods of 
radii R r j the outer radii of the core is denoted by R P j and the distance to the nearest neighboring 
rod is Rhj- for the coated rod assemblage it is clear that Rhj is determined by R r j together with 
9r and 6p. 




Figure 3: Doubly Coated Cylinder 




Figure 4: Cross section of the unit period filled with coated rods. The coated rods fill all space. 



Me// 




Figure 5: the relation between /x e // and (u/q/cY 



This type of configuration is known in the composites literature and is referred to as a doubly 
coated cylinder assemblage [35], [37] , [35], [3DJ. The explicit formula for fi e ff follows from (12.32)) 
and is given by 






(3.1) 



j=o 



m— 



jm 



The poles of \x B f / (ujq / c 2 ) are given by fc 2 m , see Figure [5] 

The explicit formula for the effective dielectric permittivity follows from a standard calculation 
using (|2~^T) . ([231]) , (13^3j) and (|3^31) (see, for example [3D]) and is given by 



i,-i, 2(1 - 9 H )€£ 



€ eff ~ 6 H 



2 ^ 



(3.2) 



with 



-l _ ,-i i 20R e p 

2&R — @P 



(3.3) 



ll<? 



On writing e p (u;g/c ) = 1 §j^ and en = 1 we obtain the formula for the frequency dependent 

effective dielectric tensor given by 

„ . / l+g«+gp I (l (^p/c) 2 ^ 1 ( 6 P \ 
l-{0R + 0p) l-(0R+6p) ^ y 1 (ua/c) 2 ) \26 n +e P ) 



eff \c 2 )~ i + e n + 9 P I i-(e n +g P ) i a t^/c)^ " 1 / g£_ \ r 

l+0fl+9p + 1/ (^,/c) 2 ) \28 R +e P/ 

The dependence of // e // and e e // on the rod diameter and coating thickness and host area fraction 
is explicitly given by (|3.ip and (|3.4p . This explicit dependence allows for the tuning of the dispersion 
relation 



where 



fc 2 = ^nt ff , (3.5) 

For each coated rod with radius R r j, outer coating radius R P j and nearest neighbor distance Rhj 
the ho field inside the rod is given by (|2.31[) and the hi field outside the rod is 

hi = (1 + R y l I rcos7 + — ^ cos 7 I , R rj < r < R pjl (3.7) 



20 R + 9p + 9 P [l 

-rcos7 + 



2 / 2 \ — 1 



F 

26 R + 9 P -9 P (l-^ I ) R 2 
+ ^ ^cj ^ cog ^ ^ < r < ^ (3 g) 

where (r, 7) are local polar coordinates inside the coated rod and F is given by 

F={26 R + 6p){l + 6 R + 6p)+[\-^i^\ 6p{\-{e R + 6 P )). (3.9) 



It is clear that ho diverges inside the rods with radii R r j as Wg/c 2 approaches k 2 m while hi diverges in 
the neighborhoods of all the coated rods, R P j < r < Rhj as Wg/c 2 approaches the zero of e e // ( "jr ) ■ 

4 Tunable double negative behavior 

When e e ff and A* e // have the same sign it is clear that n 2 e , * > and wave propagation occurrs. In 
this section, we show that it is possible to tune the assemblage through the choice of rod radii, coating 
thickness and host volume fraction to exhibit single negative, double negative, and double positive 
behavior depending on the signs of e e // and /i e //- For coated cylinder assemblages it is evident 
that the distance between neighboring rods and the coating thickness is explicitly determined by the 
area fraction of the host phase, the area fraction of the rods and the rod radii R r j . The dielectric 
function (|3.4|) has only one pole and one zero and these depend explicitly on the geometry of the 
coated rod assemblage and are given by 

D-f^ya 9 H {29 R + 8 P ) 

p y c ' e H {2d R + ep) + {2-e H )0p v ' ; 



and 



p . = f gpy» (2-9 H )(2e R + 9 P ) 

p K c> (2-9 H )(29 R + e P ) + e H e P - y ' ; 



e fff 




(^) 2 



Figure 6: the relation between e e // and (u>o/cy 



Observe that p < p* and e e // is negative when p < (ujq/c) 2 < p* . The relation between e e // and 
(loq/c) 2 is displayed in Figure [5J 

From Figure [5]and Figure [Sjit is clear that both e e ff and /i e // are negative for p < (u>q/c) 2 < p* 
provided that (ujq/c) 2 is simultaneously greater than but close to the poles of /i e //- Therefore for 
e e ff and A* e // to be simultaneously negative it is required that p < k 2 m < p*, i.e., the assemblage 
must contain rods with radii R r j satisfying the condition 



Xor, 



(2 - d h )(2d r + 9 p ) + d h t 



e r ui 2 



(2-9 h )(29 r + < 



< R r j < 



Xor, 



9 h (29 r 



e r u 2 



6h(26 r + 9 P ) 



\)9 p 



(4.3) 



Now we fix the area fractions Or = 0.8, 9p = 0.1, 9h = 0.1 , and set 



25.8 THz and 



e r = 200m 2 so that p = 0.349 x 10 1 
k2 = (7) 2 ^// 



and p* = 0.737 x 10 1 



The dispersion relation 



is displayed in Figure [7] for an assemblage with the two largest rod radii chosen to be 
R r o — 0.2 and R r \ = 0.17. These radii satisfy the requirement given by (14.31) while the radii of all 
other rods in the assemblage are chosen to lie below these two values and do not satisfy (J4.3I) . It is 
clear from the figure that, e e // < when p < (luq/c) 2 < p* and otherwise e e // > . In this example, 
p is below fc 2 and p* lies between fcg and k 2 . Notice in Figure [7] that the slope is negative for 
the dispersion curve lying between k^ and p* , indicating a negative group velocity. 



k 2 






/A 

"ill! 
pVV 



(wo/c 



Figure 7: Dispersion relation k 2 = (^a'> 2 '" 2 

THz. The assemblage contains rods of radii R t q — 0.2 and R r \ 



n^ff for a material with e p 



0.17. 



<R 



where w. 



7' 



25.8 



In Figure [51 we choose the largest two radii to be given by R r $ = 0.2 and R r \ = 0.18 and all 
other rods have radii below these values. For this case, p is below /cq but p* is between k 2 and k 2 -^. 

Notice that for points x inside the coating phase and for (wq/c) 2 close to k 2 m , that (|2.31l) and 
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(wo/c 



k 2 

Kg---- 



/I 



«o6 



^7? 



where w = 25.8 



Figure 8: Dispersion relation k 2 — (^) 2 n 2 e ff for a material with e p = 1 

THz. The largest two rod radii appearing in the assemblage are given by R r o = 0.2 and R r \ = 0.18. 



(I2T43)) show that if(x,i) behaves like h(x)e i( - kk ' x ' tu) where 

(w /c) 2 - k 2 m 



h(x) = Me " f >e Jtex (l - 



'l|V'0m))V , 0m- 



(4.4) 



In particular, when (cjo/c) 2 is near fc 2 , the profile of h is V'oo (see Figure [H] with R t q = 0.2 ) and 
when (wq/c) 2 is near fcgi , the profile of /i is V'oi (see Figure [TU]with R r Q = 0.2 ). 




Figure 9: the graph of ipi 



on 




Figure 10: the graph of ijjqi 
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